Equation of state dependence of Mach cone like structures in Au+Au collisions 
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In jet quenching, a hard QCD parton, before fragmenting into a jet of hadrons, deposits a fraction 
of its energy in the medium. As the parton moves nearly with speed of light, much greater that the 
speed of sound of the medium, quenching jet can generate Mach shock wave. We have examined the 
possibility of Mach shock wave formation due to jet quenching. Assuming that the deposited energy 
quickly thermalize, we simulate the hydrodynamic evolution of the QGP fluid with a quenching 
jet and subsequent particle production. Angular distribution of pious, averaged over all the jet 
trajectories, resembles 'conical flow' due to Mach shock wave formation. However, speed of sound 
dependence of the simulated Mach angles are at variance with that due to shock wave formation in 
a static medium or in a medium with finite velocity. 

PACS numbers: PACS numbers: 25.75.-q, 13.85.Hd, 13.87.-a 



I. INTRODUCTION 

From the general theoretical grounds, it was predicted 
[l| that in a dense, deconfined medium, high-speed par- 
tons will suffer energy loss, significantly modifying the 
fragmentation function, which in turn will lead to sup- 
pressed production of hadrons. The phenomena called 
jet quenching, was later verified at Relativistic Heavy 
Ion Collider (RHIC), in Au-|-Au collisions at ^/spfj^=200 
GeV 0-3 • While jet quenching is verified, how the lost 
energy is transported in the dense medium is uncertain. 
It has been suggested that a fraction of lost energy will 
go to collective excitation, call the "conical flow" 0-0] ■ 
The parton moves with speed of light {cjet ~ 1), greater 
than the speed of sound of the medium (c^), and the 
quenching jet can produce a shock wave with Mach cone 
angle, 9m = cos^^^p^. Resulting conical flow will have 
characteristic peaks at (j) — t: — 9m and </> = t: + 9m. Both 
in STAR Q and PHENIX [9] experiments, indication of 
such peaks are seen in azimuthal distribution of secon- 
daries associated with high px trigger in central Au-|-Au 
collisions. As Mach cone is sensitive to the speed of sound 
of the medium, it raises the possibility of measuring the 
speed of sound of the deconfined matter of Quark-Gluon- 
Plasma. Mach like structure (splitting of away side peak) 
can also be obtained in various other models, e.g. gluon 
Cerenkov like radiation models [Tol. the parton cas- 
cade model [13] ■ Recently in |l3| energy density wake 
produced by a heavy quark moving through a strongly 
coupled N=4 supersymmetric, Yang-Mills plasma is com- 
puted using ADS / CFT correspondence. Mach like struc- 
tures is also observed for quark velocity greater than the 
speed of sound of the medium. 

Mach like structures are best explored in hydrodynam- 
ics models. Several authors [14| - |25| have explored the 
possibility of Mach shock wave formation in a jet event. 



In order to understand the dynamics of energy trans- 
port by the QGP fiuid, in ^21r.-25j . we have modeled the 
quenching jet as a (time dependent) source and solved hy- 
drodynamical energy-momentum conservation equations. 
It is implicitly assumed that the energy deposited lo- 
cally along the trajectory of the jet is quickly thermal- 
ized and evolve hydrodynamically. Explicit simulations 
of the "'Hydro-I-Jet'" model [11], indicate that in a single 
jet event, azimuthal distribution of pions do not resemble 
the structure characteristic of 'conical flow'. Instead of 
two peak distribution as expected in a conical flow, the 
distribution exhibits a single peak structure. The peak 
position depend on the jet trajectory. However, when 
averaged over all the jet trajectories (di-jet can be pro- 
duced anywhere on the surface of the fireball and since all 
positions are equally fikely, one needs to average over all 
possible di-jet production points) , hydrodynamical evolu- 
tion with a quenching jet produces a angular distribution 
characteristic of 'conical flow, i.e. two peaks separated by 
a dip at (j) = TT [2^ . When simulated with realistic initial 
condition, equation of state etc, azimuthal distribution 
of excess pions, averaged over all possible jet trajecto- 
ries, reasonably well reproduce the STAR and PHENIX 
data [^ on the dijet-hadron correlation in 0-5% central- 
ity Au+Au collisions [22l - [25| . However, it is yet uncer- 
tain, whether or not the observed dijet-hadron correla- 
tion is due to Mach shock wave formation. In an evolving 
medium, due to flnite fluid velocity, Mach shock fronts 
are distorted [2Q]. Model calculations [1^ indicate that 
the inside shock front is pushed in and the outside shock 
front is pushed out. If 9^.j and 9'^^ denote the outside 
and inside shock cone angles, then for fiuid velocity u, 
the Mach angles are changed as [1^ , 
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where 7u = (1 — u^) 7s = (1 — c^) For super- 
sonic flow i.e. at u > Cs, the inside cone angle 9 
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and jet trajectory lies outside the Mach region. Then de- 
pending upon the fluid velocity, inside shock front may 
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FIG. 1: schematic representation of a jet moving through the 
medium. The high pr pair is assumed to produce at P on the 
surface of the fireball characterized by the angle (j)prod- One 
of the jet escapes forming the trigger jet, the other move in 
the fireball at an angle (pjet- 



the fireball. For Au+Au collisions at impact parameter 
b, the fireball boundary is an ellipsoid with minor and 
major axis, A = R - b/2 and B ^ R^l - 62/4i?^ with 
R=6.4 fm. As shown in Fig|TJ the jet production point 
can be characterized by the angle ipprod (0 < (j^prod < 27r) 
only . One of jet moves outward and escapes, forming 
the trigger jet. The other enters into the fireball. The 
trajectory of the jet can be designated by the angle 4>jet, 
(— ^ < (pjet < +f )• The fireball is expanding and cool- 
ing. The ingoing parton travels at the speed of light and 
loses energy in the fireball. We note that just at the cre- 
ation point P, the jet is a hard "undressed" parton. The 
"undressed" parton loses energy by constantly emitting 
radiation (gluons) which in turn emit new radiation, de- 
veloping a shower. The shower has a complicated struc- 
ture, but multiplicity in the shower grows non-linearly 
and very soon the region near the head of the jet can 
become a "macroscopic" body providing large perturba- 
tion to the medium. We treat this perturbation with 
a simple assumption. We assume that the perturbation 
is thermalized and acts as a source of energy and mo- 
mentum to the fireball. We solve the energy-momentum 
conservation equation. 



be obliterated and only the outside shock front survives. 
Angular distribution due to distorted shock fronts may 
lead to single peaked angular distribution, as obtained in 
simulations in the 'Hydro-|-Jet' model [23]. On averag- 
ing over different jet trajectories, two peak structures can 
still be obtained. However, in that case, peak positions 
will corresponds to (f> — 7r/2 ± O'j^j and simulated Mach 
angles will show definite dependence on speed of sound 
of the medium, decreasing with the speed of sound. 

In the present paper, in the " 'Hydro-I- Jet" ' model [2l| . 
we have investigated the speed of sound dependence of 
the azimuthal distribution of secondaries. For a simple 
equation of state, p = c^e, azimuthal distribution of 7r~ 
is studied as a function of speed of sound, Cg=0.1-0.4. 
It is shown that the equation of state dependence of the 
azimuthal distribution of 7r~ is unlike that due to a Mach 
shock wave in a static or moving medium. 

The paper is organized as follows: in section |IT1 we 
describe the model used to study jet quenching in an 
expanding medium. In section Hill simulation results are 
discussed. Finally, summary and conclusions are given 
in section ITVl 



II. HYDRO+JET MODEL 

A schematic representation of the jet moving through 
the medium is shown in FigH] We assume that just be- 
fore hydrodynamics become applicable, a pair of high-p^ 
partons is produced. Strong jet quenching and survival 
of the trigger jet, forbids production in the interior of the 
fireball. Jet pairs can be produced only on a thin shell on 
the surface of the fireball. Without any loss of generality, 
we assume that the pair is produced on the boundary of 



where the source is modeled as. 
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J^ix) = J{x) {l, -COs{(j)jet), -sin{(j)jet),0) , (3) 
'''' ' ^^(r-r,et(t)). (4) 
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Massless partons have light-like 4-momentum, so the cur- 
rent J"^ describing the 4-momentum lost and deposited 
in the medium by the fast parton is taken to be light-like, 
too. rjot(t) is the trajectory of the jet moving with speed 
\dxjot/dt\ —c. ^{x) is the energy loss rate of the parton 
as it moves through the liquid. It depends on the fluid's 
local rest frame particle density and we parameterise it 
as. 
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where s{x) is the local entropy density without the jet. 
The measured suppression of high-p^ particle production 
in Au-I-Au collisions at RHIC was shown to be consis- 
tent with a parton energy loss of =14GeV/fm at 
a reference entropy density of sq — 140 fm~^ [23|. In [2^ . 
it was shown that with realistic initial conditions and 
equation of state, hydrodynamic evolution with parton 
energy loss of — 14GeV/fm at a reference entropy 
density of so = 140fm~^ reasonably explain the dijet- 
hadron correlation. Presently, we are interested in the 
equation of state dependence on Mach cone like surfaces. 
We have considered simple equation of state p = c^e, 
with €3=0.1,0.2,0.3 and 0.4. For the present study. 
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with a view to enhance the effect of jet quenching, we 
use, ^ Ip =28 GeV/fm at a reference entropy density of 

so = 140fm^^. Reference energy loss =28GeV/fm 
may seem to high, but energy deposited to the medium 
is never too high due to the factor s{x)/{so — 140) in 
Eqini The di-jet is produced at the peripheri where lo- 
cal entropy density s{x) << sq and at early time the 
energy loss is not large. By the time, jet has reached 
the centre, local entropy density is degraded very much. 
For example, in boost-invariant longitudinal motion ( in 
Bjorken dynamics), for a jet moving along the diameter, 
by the time jet reaches the centre, local entropy density 
is degraded by a factor of 10 and again the jet energy 
loss is not large. Indeed, we have checked that average 
energy loss is ~ 10-15 GeV, depending on the jet tra- 
jectory. It may be mentioned that even though in eqlSj 
partonic energy loss is assumed to depend only on the 
fluid rest frame density, it could as well depend on the 
parton energy [28l - [30t . 

For the hydrodynamic evolution we use a modified ver- 
sion of the publicly available hydrodynamic code AZHY- 
DRO [Slilsll. The code is formulated in (r, a;,y, 77) co- 
ordinates, where T^\/t^—z'^ is the longitudinal proper 
time, 77=iln |^ is space-time rapidity, and = (x, y) 

defines the plane transverse to the beam direction z. 
AZHYDRO employs longitudinal boost invariance along 
z but this is violated by the source term (j4]). We there- 
fore modify the latter by replacing the (5-function in (|4]) 

by 



'5'(r-rjet(t)) 



- 6{x - Xjctir)) S{y - yjct(T)) 



(6) 



with G = 0.35 fm. 

Intuitively, this replaces the 'needle' (jet) pushing 
through the medium at one point by a 'knife' cutting 
the medium along its entire length along the beam di- 
rection. If we disregard the radial motion, the motion 
of a pointlike source produces a 'conical' wave but mo- 
tion of a 'linelike' source will produce a 'wedge like' flow. 
In a 'wedge like' flow, effect of jet quenching will be en- 
hanced, the two peak angular structure will be more pro- 
nounced, than in a 'conical' flow. However, the peak 
positions themselves (or the cone angle) will remain un- 
changed. In the present study, we are interested in the 
equation of state dependence of Mach cone like surfaces 
produced by a quenching jet. While a complete study 
of jet quenching do require a full (3-t-l)-dimensional hy- 
drodynamic calculation, qualitative change in Mach cone 
like surfaces due to different equation of states can be ad- 
equately addressed in a boost-invariant simulation. We 
have tried to estimate the enahancement expected due 
to the assumption of boost invariance. PHOBOS collab- 
oration has measured the pseudo-rapidity density for the 
charged particles [33*1 . Over a pseudo-rapidity range ^ 
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FIG. 2: Contour plot of energy deposition vax — r plane, by a 
quenching jet in medium with square speed of sound Cs=0.2 
(panels (a) and (b)) and Cg=0.4 (panels (c) and (d)). The 
quenching jet is produced at (x=4.5 fm, y=5 fm) and moving 
parallel to the X-axis. Energy deposition at fixed y=4.5 fm 
and y=-3 fm are shown. Note that away from the trajectory, 
the quenching jet deposit negligible energy in to the medium. 



(-6 to 6) boost-invarinace is approximately valid in the 
range ~ (-2.5 to -|- 2.5). The assumption boost-invariance 
then overestimate the particle yield by a factor ^1.5. Ef- 
fect of jet quenching will be overestimated by a similar 
factor. In less central collisions, the effect will be still 
less overestiamted. It may also be mentioned here that 
the Gaussian like energy loss distribution (see Eq|B]) is at 
variance with realistic calculations. Microscopic studies 
[33 - [36| do indicate that the parton energy loss is more 
like a 'spike' at short distances and an exponential fall 
off at large distances. 

The modified hydrodynamic equations in (r, a;, y, 77) co- 
ordinates read |2l| 

d,f^^ ^d^(vj-'^)^dy{vyf^^) = ~d^v~h (8) 

drf^y ^ d^Av^^f^y) + dy{Vyf^y) = -dyP, (9) 

where f^'•'^TT^"', -D, =T"/r^^, p^rp, and J = tJ. 

To simulate central Au-|-Au collisions at RHIC, we use 
the standard initialization described in 13 ill and provided 
in the downloaded AZHYDRO input file [32| , correspond- 
ing to a peak initial energy density of Si = 30 GeV/ fm^ 
at Ti = 0.6 fm/c. Since we are exploring the equation of 
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FIG. 3: (color online) The black lines are transverse momen- 
tum spectrum of n~ from evolution with out a quenching jet. 
The blue lines are from evolution with a quenching jet. 

state dependence on the splitting of the away side jet, 
we assume a simple equation of state, p{e) = c^e, with 
e{T) = ^gqT'^, gq = 47.5. We also neglect the possibility 
of phase transition. Till the freeze-out the medium evolve 
as QGP. We assume instantaneous hadronisation at the 
freeze-out temperature Tp=100 MeV and using Cooper- 
Frye prescription compute invariant distribution of tt^. 
It must be mentioned that in reality, QGP fluid will cool 
with time and at a critical temperature (Tc) will undergo 
a phase transition to hadronic matter. However, speed of 
sound dependence of Mach cone will then be more com- 
plicated. If Ti and Tp are the initial and freeze-out time 
then in an evolving medium the Mach cone angle will be 
defined as 

1 r'" 

COsOm = / Cs{T)dT (10) 

If Cs change between and Tp, unraveling speed of 
sound dependence of the Mach like structures will be 
more complicated. 

III. RESULTS 

In the following, we will simulate Au+Au collisions at 
zero impact parameter. In b=0 Au+Au collision, the re- 
action zone is spherically symmetric. For a spherically 
symmetric fireball, the fireball can be rotated by the an- 
gle (pjet to make <j)jet=^ and the jet trajectory can be 
characterized by the angle 4>prod only, (jyprod varying be- 
tween [— 7r/2, 7r/2]. 

Let us discuss speed of sound dependence of energy de- 
position by a quenching jet. In Figl2l we have shown the 
contour plot of energy deposition by a quenching jet in 
medium with square speed of sound c^=0.2 and 0.4. The 
quenching jet is produced at (x=4.5 fm, y=5 fm) and is 



moving inward parallel to the x-axis, and the contours 
are drawn in the x — t plane at a fixed y=5 fm (panels 
(a) and (c)) and y=-3 fm (panels (b) and (d)). With 
time, the quenching jet deposits less and less energy. It 
is understood. Jet energy loss is weighted by the entropy 
density. As the fiuid evolves, entropy density drops and 
the jet deposit less and less energy to the medium. One 
also notice that away from the jet trajectory, the jet de- 
posit hardly any energy in the medium. For example, 
energy deposition falls by a factor of ~ 10~^^ from y=5 
fm to y=-3 fm. It is interesting to note that the perturba- 
tion to the medium with cl—Q.2 or 0.4, by the quenching 
jet, is approximately same. The quenching jet deposit 
approximately same energy to the medium. Other con- 
ditions remaining the same, low medium evolve for 
longer duration than in high Cs medium. However, the 
energy transfer from jet to medium is inefficient at late 
time and approximately same energy is deposited in the 
medium. 

It is instructive to compare transverse momentum dis- 
tribution of secondaries from evolution of fiuid with and 
with out a quenching jet. The black solid, long dashed, 
medium dashed and short dashed lines in Figl^l are the 
transverse momentum distribution of tt^ , from fluid with 
square speed of sound, c^^O.l, 0.2, 0.3 and 0.4 respec- 
tively. Other conditions remaining the same, particle 
production depend considerably on the speed of sound. 
With increasing speed of sound, px spectra flattens, but 
yield at low px decreases. Freeze-out occur early in high 
Cs medium than in a low Cs medium. Flattening of the 
px spectra is indicative of early freeze-out. If px spec- 
tra is approximated as dN/cPpx oc exp{—Mx/T), effec- 
tive source temperature increases 'linearly' with speed of 
sound. In FiglSj the blue lines are transverse momentum 
distribution of tt" with a quenching jet. The jet pro- 
duction angle is 0prod=l-O5 rad. Effect of jet quenching 
is not large on the transverse momentum distribution. 
Yield at high px is increased by 10-20%. At low px, 
effect of the quenching jet is even less. 

Effect of jet quenching is better observed in azimuthal 
distribution of n~ . In Fig|31 azimuthal distribution of 
7r~ (0 < Pt < IGeV) from evolution of fluid with and 
with out any quenching jet are shown. The black solid, 
long dashed, medium dashed and short dashed lines are 
azimuthal distribution of 7r~ form evolution of fluid with 
square speed of sound 0^=0.1, 0.2, 0.3 and 0.4 respec- 
tively. The blue lines are same from fluid evolution with 
a quenching jet. Without any quenching jet, the pion 
distribution is flat. In b=0 Au-|-Au collisions, reaction 
zone is azimuthally symmetric and so is the azimuthal 
distribution of tt". With a quenching jet, tt" distribu- 
tion is modified. The quenching jet defines a direction 
and azimuthal symmetry is lost. Pion distribution is 
no longer flat but shows a single peak structure. The 
peak positions appear to depend marginally on the speed 
of sound. We also note that in fluid evolution with a 
quenching jet, 7r~ production is depleted around w 1 
rad. In FigUl the arrow marked the jet production angle 
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FIG. 4: (color online) The black lines are azimuthal distribu- 
tion pions from evolution of QGP fluid with square speed of 
sound c^=0. 1,0.2,0.3 and 0.4. The blue lines are same with a 
quenching jet. The jet is produced at 0prod=l.O5 rad, indi- 
cated b ' ' 
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FIG. 5; same as in Fig|4]but in the pt range IGeV < pt < 
2GeV. 



4>prod- 7r~ production is depleted approximately at the 
entry point of the jet. The azimuthal distribution of 7r~ 
with a quenching jet is unlike that expected in conical 
flow due to shock wave. In conical flow, azimuthal dis- 
tribution is expected to be double peaked with a dip at 
(j) — TT. However, as seen here, distribution is hardly mod- 
ified at (f) > TT. Also, variation of peak angle with Cg is 
much less and opposite to that expected in a Mach shock 
wave. In Fig (5] and [6l azimuthal distribution of tt^ in px 
range, IGeV < pr < 2GeV and 2GeV < pr < 3GeV are 
shown. Qualitatively, the azimuthal distribution of high 
Pt pions are similar to that of low px pions. At high 
Pt also, azimuthal distribution is single peaked, associ- 
ated with depleted production around the jet production 
angle. However, the peaks are better defined at high px- 
In Fig(71 we have mapped the azimuthal distribution of 



FIG. 6: same as in Fig[4]but in the pr range 2GeV < pt < 
3GeV. . 



""^ (1 < Pt < 2 GeV) in x-y plane. Remember that in a 
dynamic model like hydrodynamics, the particles are con- 
tinuously emitted from a dynamic freeze-out surface. The 
mapping x = ^cos{(j)), y = ^sin{4>) give the shape of 
the effective freeze-out surface. Black lines in Figl7]corre- 
spond to fluid evolution without any quenching jet. The 
freeze-out surface is a circle then. Effective freeze-out 
surface is modified with a quenching jet. The modifica- 
tion is qualitatively similar for all values of c^. The sur- 
face, in the lower half plane is largely unaltered. Surface 
on the upper half plane is modified; around the jet entry 
point, the surface is pushed in. As a consequence, the 
surface is pushed out in another direction. Apparently 
the matter behaves as an incompressible fluid. Modifi- 
cation of surface is also great when speed of sound of 
the medium is large. However, the modification is un- 
like that expected in a Mach shock wave formation. The 
simmulation studies do indicate that in a single jet event, 
quenching jet do not lead to a double peak structure in 
azimuthal correlation. The result is in agreement wit re- 
cent studies of azimuthal correlation induced by energy 
loss by a heavy quark [l^. It was shown that a pQCD 
source term for heavy quark energy loss do not lead to 
double peak structure in azimuthal correlation. 

Even though, in a single jet trajectory, azimuthal dis- 
tribution of TT^ do not show a double peaked structure as 
in a conical flow, when averaged over the jet trajectories, 
a double peaked structure, mimicking the conical flow, 
can emerge ^22,, 2^ ■ A di-jet can be produced anywhere 
on the surface of the fireball and one need to average 
over all the jet trajectories. As demonstrated in (2^, a 
quenching jet in the upper half plane produces excess 
pions in the angular range < (/) < tt, while a jet in 
the lower half plane produces excess pions in the angular 
range tt < < 27r and on averaging produces a double 
peaked azimuthal distribution with a dip at = tt. 

In FigjSl we have shown the azimuthal distribution 
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FIG. 7: (color online) Azimuthal distribution of 7r~ is mapped 
into x-y plane. The black and red lines in panel (a), (b), (c) 
and (d) corresponds to fluid evolution without any quenching 
jet and with a quenching jet produced at 0=0.78 rad. 




of 'jet trajectory averaged' excess pions [y^ . ^ 
^ . Trajectory averaged tt^ distribution is ob- 
tained as, 



dN 



jet 



1 dN{cj)prod) ^ 
-7r/2 d(j) 



i^prod 



(11) 



Excess pion distribution in the px range, (a) < Pt ^ 
IGeV, (b) IGeV < pr < 2GeV and (c) 2GeV < pr < 
3GeV are shown in three panels of FiglH Jet trajectory 
averaged azimuthal distribution clearly show a two peak 
structure with a dip at = tt. At high px and also in 
medium with high speed of sound, the peak structures 
get better defined and the dip at (p—ir deepens. We also 
notice that in the simulation, due to quenching jet, par- 
ticle production is depleted around (j)—0, i.e. around the 
trigger jet. Depletion around trigger jet is not observed 
experimentally. We have neglected the trigger jet. The 
trigger jet fragment in the vacuum and populate phase 
space around (f> = 0. Particles from evolution of the 
fireball and particles from trigger jet fragmentation can 
not be distinguished and depletion of particle production 
around = is not observed. 

Azimuthal distribution of trajectory averaged excess 
pion distribution resemble a conical flow, i.e. two peak 
distribution with a dip at (j) = tt. As argued earlier, if 
the peaks are due to distorted shock waves in a mov- 
ing medium, such that inside shock front is obliterated 
and only the outside shock front survives, then the az- 
imuthal distribution of excess pions with have peaks at 
4>peak = TT ± 6*^. Peak position or the outside shock cone 
angle O'^j will show definite variation with speed of sound 
(see EqlT]). To test this, we have fitted the azimuthal dis- 
tribution of excess pions by a double Gaussian, 



FIG. 8: Jet trajectory averaged azimuthal distribution of ex- 
cess pions [^liet - ^Inojet] , iu the pT range, (a) < pt < 
1 GeV, (b) IGeV < pr < 2GeV and (c) 2GeV < pr < 3GeV. 
The blank circles, squares, up triangles and down triangles 
corresponds to square speed of sound 0^=0.1, 0.2, 0.3 and 
0.4 respectively. The lines are fit to the simulated azimuthal 
distribution by two Gaussians (EafT2]). 



dN 



cx e 



(12) 



Fits obtained to the simulated tt^ distribution is shown 
in FigEl Double Gaussian fits the azimuthal distribu- 
tion of the trajectory averaged excess pions due to jet 
quenching. In Fig|9l we have compared the simulated 
Mach angles (O'j^/^"^) with Mach angles (Eq[T]) in a mov- 
ing medium. The solid lines in Fig|9]are Mach angles in a 
moving medium, as a function of square speed of sound. 
We have shown the angles for fluid velocity, tt=0, 0.2, 
0.4 and 0.6. Mach angles in a moving medium decreases 
with speed of sound and also with increasing fluid veloc- 
ity. The filled circles, squares and triangles in FiglS] are 
simulated Mach angles 9^/"" from the excess pion distri- 
butions in the pT range < < IGeV, 1 < pt ^ 2GeV 
and 2 < pt < 3GeV respectively. Unlike in theory, sim- 
ulated Mach angles depend on pt- Speed of sound de- 
pendence of d'l/^"^ is also inconsistent with the speed of 
sound dependence of Mach angles in a moving medium. 
Simulated Mach angles will be consistent with Mach an- 
gles in a moving medium, if, on the average, fluid velocity 
is less in medium with high speed of sound. For exam- 
ple, if fluid velocity in medium with 0^=0.1, 0.2, 0.3 and 
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FIG. 9: The solid lines (from top to bottom) are theoretical 
estimate of Mach angle [see Eq[T], as a function of the 
square speed of sound, for fluid velocity u=0, 0.2, 0.4 and 
0.6 respectively. Simulated Mach angles (Sj^^™) from jet 
trajectory averaged vr" distribution, in the pr range, < 
PT < 1 GeV, 1 < PT < 2 GeV and 2 < pr < 3 GeV are 
shown as filled circles, squares and triangles. 

0.4 are ^ 0.5, 0.38, 0.25 and 0.1 respectively, simulated 
Mach angles in the px range < pt < IGeV, will be 
consistent with the Mach angles in a moving medium. 
However, in hydrodynamic evolution, fluid velocity grow 
faster in a high Cs medium than in a low Cg medium. 
Consistency condition i.e. fluid velocity is less in high 
Cs medium than in low Cs medium can not be satisfied. 
Speed of sound dependence of the simulated Mach angles 
in the '"Hydro+Jet"' model, can not be reconciled with 



theoretical estimate of Mach angles in a fluid with finite 
velocity. The result agree with the finding in 18} . In [ij] 
it was shown that supersonic strings in AdS/CFT cor- 
respondence, do not produce observable conical angular 
correlation in strict Nc — >■ oo supergravity limit. How- 
ever, a special non-equilibrium neck zone, near the jet, 
do produce Mach cone like peaks, which does not follow 
Mach's law. 



IV. SUMMARY AND CONCLUSIONS 

Possibility of Mach shock wave formation due to a 
quenching jet is studied function of the speed of 

sound of the medium. Modeling the quenching jet as a 
source, we have solved energy-momentum conservation 
equation of QGP fluid with a simple equation of state, 
p — c^e, for Cg=0.1, 0.2, 0.3 and 0.4. Phase transition is 
neglected. It is assumed that at the freeze-out tempera- 
ture rF=100 MeV, QGP fluid instantly hadronizes. In a 
single jet event, simulated azimuthal distribution of pions 
due to a quenching jet does not resemble a conical flow. 
In a conical fiow, a double peak structure, separated by 
twice the Mach angle {9m — cos^^Cs) is expected. Unlike 
in a conical fiow, in a single jet event, simulated distribu- 
tion is single peaked, the peak position varying little with 
speed of sound. When averaged over all the jet trajec- 
tories, a double peaked structure, reminiscent of 'conical 
flow' emerges. However, the peak positions are at vari- 
ance with that expected in a shock wave formation. We 
conclude that the fine structure observed in experimen- 
tal azimuthal distribution of secondaries associated with 
high pt trigger in Au-|-Au collisions are possibly unre- 
lated from Mach shock wave formation. 
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